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1. 


INTRODUCTION 


Instrument  technology  developments  and  the  desire  for  "all  weather"  capa¬ 
bilities  have  contributed  to  an  increased  interest  in  the  design  and  deploy¬ 
ment  of  microwave  and  millimeter  wave  communications  and  sensor  systems. 
Nowhere  is  this  trend  more  apparent  than  in  the  meteorological  sensor  comple¬ 
ment  of  the  current  Defense  Meteorological  Satellite  Program  (DMSP)  Block  5D-3 
which  consists  of  three  microwave  instruments,  the  SSM/T  temperature  sounder, 
the  SSM/T-2  moisture  sounder,  and  the  SSM/I  imager  (Isaacs  et  al ,  1986; 

Falcone  and  Isaacs,  1987).  The  follow-on  system  to  the  operational  TIROS 
sounding  system  the  ATOVS  (Advanced  TIROS  operational  vertical  sounder)  will 
also  have  a  microwave /millimeter  wave  component,  the  AMSU  (Advanced  microwave 
sounding  unit) .  There  is  also  increased  interest  in  ground  based  microwave 
sensor  systems  and  in  the  application  of  microwave /millimeter  wave  sensors  for 
airborne  surveillance.  These  developments  emphasize  the  need  for  an  accurate, 
readily  available  microwave  transmittance ,  attenuation,  and  brightness 
temperature  code  for  Air  Force  and  other  Department  of  Defense  users. 

The  present  RADTRAN  code  has  been  widely  applied  within  the  DoD  and  else¬ 
where  to  provide  the  capability  to  perform  sensor  simulation  studies  for 
satellite  borne  and  ground  based  microwave  and  millimeter  wave  systems.  In 
the  course  of  our  own  remote  sensing  related  research,  AER  investigators  have 
had  the  opportunity  to  successfully  apply  the  RADTRAN  model  in  a  variety  of 
contexts.  Most  recently,  for  example,  we  have  studied  the  effect  of  clouds  on 
millimeter  wave  moisture  soundings  using  a  modified  version  of  RADTRAN  for  our 
sensor  simulations  (Isaacs  and  Deblonde,  1987). 

Our  experience  has  indicated  that  RADTRAN,  in  its  present  form,  suffers 
from  a  number  of  minor  deficiencies  which  decrease  its  general  applicability 
and  user  friendliness.  To  circumvent  these  problem  areas,  we  have  developed 
specific  solutions  for  application  to  our  in-house  RADTRAN-based  calcula¬ 
tions.  The  approaches  which  we  have  adopted  are  relevant  to  applications  of 
the  general  user  and,  therefore,  we  believe  would  significantly  enhance  the 
potential  of  the  algorithm  for  general  use. 

In  this  report  we  describe  the  physical  basis  of  these  enhancements  and 
their  implementation  in  the  RADTRAN  algorithm  in  detail.  We  describe  the 
current  RADTRAN  algorithm,  enumerate  its  deficiencies,  and  our  enhancements  to 
the  RADTRAN  algorithm  to:  (a)  update  the  code  structure  providing  user 
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friendly  menus,  internal  comments,  and  ANSI  standard  FORTRAN  77  coding, 

(b)  provide  for  consistency  with  standardized  input  model  atmospheres 
currently  implemented  within  other  AFGL  transmission/radiance  models,  (c) 
undertake  a  comparison  of  the  current  algorithm  with  contemporary  state-of- 
the-art  transmittance/radiance  codes  to  insure  its  fidelity  and  representa¬ 
tiveness,  (d)  enhance  RADTRAN  to  treat  multiple  scattering,  (e)  provide  for 
internal  calculation  of  polarized  surface  emissivities  based  on  user  specified 
surface  types,  and  (f)  document  the  code.  Finally,  we  include  appropriate 
statistical  retrieval  schemes  so  that  users  can  perform  sensor 
simulation/retrieval  studies  if  desired. 

2 .  BACKGROUND 

The  RAbi'RAN  computer  code  was  developed  by  the  Air  Force  Geophysics 
Laboratory  (AFGL)  to  provide  atmospheric  attenuation  and  brightness  tempera¬ 
ture  calculations  for  typical  atmospheric  paths  over  the  frequency  range  from 
1  to  300  GHz  (Falcone  et  al.,  1982).  This  code  has  been  utilized  and  tested 
against  experimental  data  for  frequencies  up  to  1000  GHz.  RADTRAN  provides  a 
design  tool  which  can  readily  be  used  to  assess  potential  environmental 
impacts  on  microwave  and  millimeter  wave  sensors.  The  atmospheric  attenuation 
submodels  of  the  clear  atmosphere,  fog,  cloud,  and  rain  used  in  RADTRAN  have 
been  thoroughly  documented  (Falcone  et  al.,  1979).  At  the  frequencies  of 
interest,  clear  sky  absorption  is  due  primarily  to  water  vapor  and  oxygen 
(Waters,  1976).  Ozone  also  absorbs  in  this  region,  but  it  generally  has  a 
negligible  effect  on  brightness  temperature  simulations  of  defense  and 
commercial  interest.  The  absorption  of  water  vapor  in  RADTRAN  is  evaluated 
using  the  expression  of  Barret  and  Chung  (1962)  for  frequencies  less  than  60 
GHz.  For  frequencies  between  60  and  300  GHz  the  absorption  is  evaluated  from 
the  183  GHz  line  plus  the  nonresonant  background  using  the  model  of  Gaut  and 
Reifenstein  (1971).  At  higher  frequencies,  water  vapor  absorption  is  modeled 
using  the  Van  Vleck-Weisskopf  (1945)  line  shape,  a  set  of  54  rotational  lines 
(see  Table  1)  and  the  Van  Vleck  frequency  fitting  to  the  continuum.  A  more 
recent  model  is  the  Gaut-Reifenstein  continuum  which  provides  an  adequate 
model  at  moderate  relative  humidities  (Liebe  and  Layton,  1983)  but  has  ques¬ 
tionable  temperature  dependence  (Burch,  1981).  Oxygen  absorption  is  evaluated 
using  the  parameters  of  Meeks  and  Lilley  (1963).  The  effects  of  line  coupling 
at  60  GHz  (Rosenkranz,  1975)  are  included  as  an  option  in  the  current  RADTRAN 
code.  Many  of  these  expressions  are  summarized  in  Falcone  et  al .  (1971). 
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The  current  RADTRAN  code  calculates  constituent  dependent  attenuation 
(e.g.  oxygen,  water  vapor,  clouds,  rain),  total  attenuation,  transmittance, 
and  brightness  temperature  as  functions  of  frequency  and  atmospheric,  surface, 
and  spectral  input  data.  The  brightness  temperature  calculation  is  based  on 
evaluation  of  the  radiative  transfer  equation  for  thermal  emission  (i.e. 
neglecting  scattering  sources)  of  microwave  frequencies: 

p 

Tb„  “  [<STS  +  (1  -  «a>  J  *  T(p)dr;]  Mp,.)  +  J°  T(p)dr„  (1) 

O  P<! 


where 


P 

Tu( P)  “  e>«P  I-  fQ  k  (p,p')dp'/M) 


(2) 


and 


T J  (P>  -  exp  (-  Jps  k(i/,p'  )dp'/M)  - 


(3) 


Here,  Tfl  is  the  brightness  temperature,  p  is  the  cosine  of  the  path  zenith 
angle,  and  and  are  the  upward  and  downward  transmission  functions,  re¬ 
spectively.  Dependence  on  the  temperature  profile,  surface  temperature  and 
surface  emissivity  are  explicit  in  the  above  expression.  Dependence  on  rele¬ 
vant  constituent  abundances  (such  as  oxygen  and  water  vapor)  and  spectral  line 
shape  parameter*  are  implicit  in  the  evaluation  of  the  transmittance  functions 
(Equations  2  and  3)  through  the  absorption  coefficients.  These  parameters  are 
also  functions  of  temperature  and  pressure.  The  surface  emissivity,  eg  is 
chosen  for  surfaces  of  interest  or  calculated  by  the  Fresnel  equations. 

RADTRAN  also  provides  a  capability  to  include  the  attenuation  of  cloud, 
fog,  and  rain  in  the  evaluation  of  atmospheric  transmittance.  These  models 
are  described  in  Falcone  et  al.  (1979).  The  cloud  and  rain  subroutines,  GMFOG 
and  GMRAIN,  respectively,  provide  cloud  models  which  are  largely  based  on  the 
data  of  Silverman  and  Sprague  (1970).  Cloud  attenuation  is  evaluated  using 
the  Rayleigh  approximation  and  the  mass  density  of  the  selected  cloud  model. 
Complex  index  of  refraction  data  for  water  clouds  is  taken  from  Ray  (1972). 
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Table  1.  Water  Vapor  Rotational  Line  Parameters 


Transition  Frequency  Parity  Matrix  Energy  Levels  Linewidth  > 
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.5091 

.  72 

2,11-2,20 

1227.8825 

OOEE 

1.2594 

95.19 

136.15 

.09792 

.4658 

.67 

6,52-7,43 

1227.9451 

OEEO 

.253 

888.70 

931.33 

.0688 

.2682 

.45 

7,34-8,27 

1294.4328 

OEEO 

.  184 

842.51 

885.69 

.0819 

.4577 

.  55 

9,18-8,45 

1309.7213 

OEEO 

.047 

1079.20 

1122.89 

.060 

.348 

.25 

5,32-6,25 

1323.2113 

OEEO 

.3117 

508.80 

552.94 

.08313 

.4939 

.  571 

9,81-10,74 

1329.8063 

EOOE 

.173 

2010.19 

2054.55 

.0390 

.2077 

.39 

9,82-10,73 

1329.8063 

EEOO 

.173 

2010.19 

2054.55 

.0390 

.2077 

.  39 

8,17-7,44 

1342.6967 

OOEE 

.036 

882.97 

927.76 

.066 

.375 

.  30 

5,14-5,23 

1407.4483 

OEEO 

4.2239 

399.44 

446.39 

.09470 

.5123 

.  722 

10,19-9,46 

1423.3364 

OOEE 

.059 

1293.22 

1340.70 

.055 

.322 

.24 

6,33-7,26 

1435.9270 

OOEE 

.258 

661.54 

709.44 

.0830 

.4642 

.59 

3.  TECHNIC M  APPROACH 

Although  the  RADTRAN  code  already  possesses  many  of  the  attributes  de¬ 
sired  of  a  generally  applicable  attenuation/transmittance/brightness 
temperature  simulation  code  for  the  microwave  spectral  region,  it  lacks  a 
number  of  desirable  and  potentially  useful  capabilities.  These  include  the 
capability  to:  (a)  model  frequency  dependent,  polarized  surface  emissivity, 
(b)  calculate  the  scattering  properties  of  precipitation,  (c)  perform  scalar 
multiple  scattering  calculations,  (d)  calculate  polarization  dependent, 
multiple  scattered  brightness  temperatures,  and  (e)  perform  sensor  retrieval 
simulations.  Figure  1  (Isaacs  et  al.,  1985)  illustrates  the  configuration  of 
the  completed  model  including  the  enhancements  alluded  to  above .  The 
schematic  in  Figure  1  is  annotated  along  the  horizontal  axis  to  define 
functional  subelements  in  the  program  corresponding  to  the  topics  discussed 
below.  For  example,  the  treatment  of  precipitation  scattering  occurs  in  the 
segment  called  physical  models,  while  that  for  the  filter  function  appears  in 
system  response.  Treatments  of  these  aspects  of  code  enhancement  are 
discussed  below.  Additionally,  the  baseline  code  has  been  upgraded  to  conform 
to  modern  programming  practices.  Finally,  we  describe  the  addition  of 
retrieval  capability  employing  statistical  retrieval  codes  developed  at  AER. 
This  addition  to  the  microwave  transmittance/brightness  temperature  code  would 
be  extremely  useful  for  satellite  and  ground  based  temperature  and  water  vapor 
profiling . 

3 . 1  RADTRAN  Program  Attributes 

The  existing  code  is  unstructured,  for  the  most  part,  with  the  essential 
calculation  contained  in  a  large  main  program.  Subroutines  are  utilized  for 
cloud  properties  and  evaluating  gaseous  absorption  parameter izations .  The 
code  is  largely  undocumented  with  some  code  segments  not  utilized  in  the 
calculation.  The  important  elements  of  the  calculation,  i.e.,  (a)  input  data 
files  (atmospheric  profiles,  line  spectral  data,  and  supplementary  block 
data);  (b)  gaseous  absorption  models;  (c)  cloud/rain  attenuation  models; 

(d)  atmospheric  layering  approach;  (e)  layer  transmittance  calculation; 

(f)  radiative  transfer  calculation;  and  (g)  output  formatting,  are  present. 
However,  efficiency  of  the  code  and  its  understandability  can  be  improved  by 
enhancing  the  code's  modularity  following  accepted  programming  practices. 
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Thi*-  calls  for  subroutining,  the  creation  of  relevant  common  blocks,  and  the 
elimination  of  unnecessary  looping  within  the  procedure.  The  code  needs  to  be 
commented  and  provided  with  documentation  including  relevant  test  cases. 

While  the  essential  elements  of  the  transmittance  calculation  and  radia¬ 
tive  transfer  code  are  inherent  in  the  model,  the  methodology  adopted  in  some 
of  the  code  procedures  should  be  evaluated  in  light  of  recent  developments. 
Some  aspects  of  the  gaseous  absorption  treatment,  for  example,  are  twenty 
years  old  and  require  updating  or  revision.  Specific  aspects  of  the  gas 
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Figure  1.  Schematic  of  proposed  microwave  transmittance/brightness 
temperature  simulation  code. 


6 


absorption  modeling  requiring  evaluation  include  the  spectral  line  data,  line 
shape,  pressure  broadening  coefficients,  line  coupling  approach,  and  treatment 
of  the  water  vapor  continuum.  In  light  of  recent  measurements  and  understand¬ 
ing  of  the  basic  theory  involved,  this  is  not  an  unreasonable  expectation. 

3 . 2  Frequency  Dependent.  Polarized  Surface  Emissivitv  Models 

The  RADTRAN  code  is  frequently  applied  to  microwave  systems  effectiveness 
simulation  studies  for  both  airborne  and  satellite-based  sensors  viewing  the 
earth's  surface.  Such  applications  require  the  capability  to  supplement  the 
atmospheric  radiative  transfer  calculation  with  realistic  models  of  the  micro- 
wave  properties  of  the  sensor  scene  field-of -view  allowing  for  variations  in 
polarization  and  frequency  dependence  of  emissivity  due  to  changes  in  surface 
type.  Frequently,  this  information  is  provided  via  available  empirical  data 
sets.  However,  the  range  of  available  field  measurements  is  limited.  In 
order  to  facilitate  the  evaluation  of  sensor  performance  throughout  a  range  of 
pertinent  geophysical  surface  scenarios,  deterministic  models  of  microwave 
surface  emission  are  desired  which  can  be  exercised  concurrently  with  the 
models  for  atmospheric  transmission.  To  achieve  this  goal,  a  set  of  microwave 
surface  emission  models  was  developed  by  AER  to  accompany  the  RADTRAN 
atmospheric  properties  model.  In  addition  to  obvious  applications  to  micro- 
wave  background  simulations,  application  of  such  models  to  remote  sensing 
simulation  studies  can  also  provide  a  theoretical  basis  for  understanding  of 
the  sensitivity  of  microwave  sensor  data  to  variations  in  relevant  hydrologic¬ 
al  parameters  such  as  soil  moisture,  snow,  and  vegetative  moisture  content. 

The  selection  of  a  simple  surface  modeling  approach  is  made  difficult  by 
the  complexity  of  geophysical  surfaces.  For  example,  homogeneous  dielectric 
slab  models  for  both  land  and  ocean  have  been  commonly  used  to  provide  the 
required  surface  emissivity  parameters  to  initiate  brightness  temperature 
simulation  calculations.  It  is  apparent  from  an  examination  of  recent  micro- 
wave  satellite  imagery  that  such  approaches  cannot  reproduce  the  fidelity  of 
the  complex  fields  of  observed  surface  properties.  This  is  due  to  the  neglect 
of  important  physical  mechanisms  such  as  scattering  by  such  approaches,  and 
their  failure  to  treat  the  spatial  inhomogeneities  in  dielectric  properties 
due  to  the  inherent  physical  structures  of  real  surfaces.  On  the  other  hand, 
it  is  necessary  to  consider  the  computational  level  of  effort  required  to 
model  all  of  the  physics,  even  if  it  were  well  understood. 
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Related  to  the  choice  of  an  appropriate  model,  is  the  detail  of  surface 
type  characterization  desired,  i.e.  how  many  surface  types  to  treat.  In 
reality,  of  course,  there  is  a  continuum  of  geophysical  surface  types.  A 
sufficiently  general  model  can  attempt  to  simulate  some  of  the  behavior  ex¬ 
hibited  by  subsets  of  surface  types  within  this  continuum  by  choosing  appro¬ 
priate  parameterizations  of  relevant  surface  properties  and  varying  them  with¬ 
in  representative  ranges.  The  choice  of  surface  types  employed  within  the 
proposed  RADTRAN  surface  modeling  package  was  based  both  on  the  desire  to 
treat  a  comprehensive  set  of  surfaces  and,  to  some  extent,  on  the  requirements 
of  potential  model  users  with  specific  surface  related  simulation  applica¬ 
tions.  The  surface  types  selected  are:  (a)  calm  and  rough  ocean,  (b)  first 
year  (FY)  and  multiyear  (MY)  sea  ice,  (c)  wet  and  dry  snow  over  land, 

(d)  moist  soil,  (e)  vegetation,  and  (f)  land.  The  land  surface  type  provides 
a  background  for  snow,  soil,  and  vegetation  models  in  addition  to  its 
potential  role  as  a  distinct  surface  type  itself.  A  menu  will  be  provided  to 
select  from  among  the  available  surface  type  choices.  Surface  types  are 
summarized  in  Table  2. 

Both  the  calm  ocean  and  land  are  modeled  as  simple  dielectric  slabs.  The 
other  surface  types,  however,  clearly  require  a  more  sophisticated  modeling 
treatment.  As  the  following  discussion  will  indicate,  it  is  not  appropriate 
to  treat  all  of  the  surface  types  delineated  above  by  a  single  formalism. 
Therefore,  two  distinct  approaches  have  been  applied  in  the  development  of 
these  surface  emission  models:  that  based  on  wave  theory  for  random  discrete 
scatterers  and  that  based  on  radiative  transfer  theory  for  continuous  random 
media  The  former  approach  is  applied  to  modeling  the  ocean  surface,  sea  ice, 
and  snow,  while  the  latter  is  applied  to  both  soil  and  vegetation.  These 
approaches  are  also  summarized  in  Table  2. 

In  the  wave  theory  approach  for  random  discrete  scatterers,  one  or  more 
layers  is  defined  consisting  of  a  dielectric  medium  with  either  uniform  pro¬ 
perties  or  containing  a  random  distribution  of  discrete  dielectric  spheres 
with  distinct  dielectric  properties.  These  latter  inclusions  give  the  medium 
scattering  properties  which  by  appropriate  choice  of  the  background  and  inclu¬ 
sion  permittivities  can  be  tuned  to  exhibit  the  observed  behavior  of  sea  ice 
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Table  2.  Surface  Model  Types  and  Modeling  Approaches 


Model 


Surface  Type  Modeling  Approach 


1 

Calm  ocean 

Dielectric  slab 

2 

Rough  ocean 

Random  discrete 

scatterers 

3 

FY  sea  ice 

Random  discrete 

scatterers 

4 

MY  sea  ice 

Random  discrete 

scatterers 

5 

Dry  snow 

Random  discrete 

scatterers 

6 

Wet  snow 

Random  discrete 

scatterers 

7 

Land 

Dielectric  slab 

8 

Wet  soil 

Continuous  random  medium 

9 

Vegetation 

Continuous  random  medium 

and  dry  snow,  for  example.  The  radiative  transfer  approach  for  continuous 
random  media  models  the  surface  from  a  different  perspective.  Some  surfaces 
are  spatially  inhomogeneous  in  their  dielectric  properties,  yet  the 
inhomogeneities  are  not  due  to  discrete  spherical  scatterers.  The  approach 
provides  an  alternative  treatment  in  which  the  permittivity  is  varied 
continuously  throughout  the  medium.  Furthermore,  these  spatial  variations  are 
parameterized  in  such  a  manner  that  the  relative  effects  of  variations  in  the 
vertical  and  horizontal  physical  structure  of  the  medium  can  be  modeled. 

These  surface  emissivity  models  are  summarized  in  a  recent  journal  article 
(Isaacs  et  al.,  1989). 

3 • 3  Scattering  Properties  of  Precipitation 

Treatment  of  precipitation  in  brightness  temperature  simulations  require 
the  specification  of  precipitation  scattering  optical  properties  including  the 
extinction  coefficient,  single  scattering  albedo,  and  the  angular  scattering 
function.  These  data  are  not  currently  implemented  within  RADTRAN  and  must  be 
supplied  by  the  user. 

Precipitation  scattering  properties  are  generally  available  via  standard 
Mie  theory  calculations.  The  Mie  theory  formalism  requires  a  knowledge  of 
particle  size  distribution  and  index  of  refraction.  The  index  of  refraction, 
in  turn,  is  dependent  on  frequency,  phase  (i.e.  ice  or  water),  and  tempera¬ 
ture.  To  avoid  the  cumbersome  necessity  of  performing  on-line  Mie  theory 
calculations  to  support  each  possible  combination  of  these  model  variables 
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within  multiple  scenario  brightness  temperature  simulations,  a  parameteriza¬ 
tion  has  been  developed  based  on  the  existing  Mie  theory  calculations  of 
Savage.  This  parameterization  is  available  for  implementation  within 
RADTRAN.  The  attributes  of  the  precipitation  property  modeling  subroutine  are 
described  in  a  paper  by  Isaacs  et  al.  (1989). 

The  resultant  subroutine  provides  an  efficient  method  to  obtain  the 
extinction  coefficient,  single  scattering  albedo,  and  angular  scattering 
function  over  the  frequency  domain  from  19  to  240  GHz.  The  angular  scattering 
function  is  given  in  terms  of  its  first  eight  Legendre  polynomial  expansion 
coefficients.  Precipitation  angular  scattering  functions  are  not  highly  an¬ 
isotropic  at  microwave  frequencies  and  this  number  of  terms  usually  suffices 
to  describe  them.  Furthermore,  this  number  of  terms  is  consistent  with  the 
Gaussian  quadrature  required  to  specify  the  brightness  temperature  field 
(Savage,  1978).  The  scattering  function  asymmetry  factor  used  in  standard 
multiple  scattering  approximations  is  easily  obtained  from  a  knowledge  of  the 
second  Legendre  coefficient. 

3.4  Multiple  Scattering 

When  precipitation  is  present  evaluation  of  the  multiply  scattered 
brightness  temperature  field  requires  a  generalization  of  the  radiance  (i.e. 
brightness  temperature)  solution.  To  generalize  the  radiance  solution,  either 
scalar  (i.e.  ignoring  polarization)  or  vector  (i.e.  including  polarization) 
multiple  scattering  treatments  are  possible.  The  scalar  approach  is  appro¬ 
priate  for  nadir  viewing  only  and  when  the  atmospheric  path  and  transmittance 
are  such  that  there  are  no  surface  contributions  to  the  brightness  temperature 
field.  (Surface  emission  and  reflection  can  be  highly  polarized  depending  on 
the  surface  type,  viewing  angle,  and  frequency.)  We  have  implemented  a  scalar 
multiple  scattering  approximation  into  FASCODE  based  on  an  approximation  of 
the  multiple  scattering  source  function  (Isaacs  et  al.,  1987c). 

A  more  complete  vector  treatment  is  required  in  the  microwave  region, 
capable  of  treating  the  multiple  scattering  of  microwaves  by  anisotropic,  in¬ 
homogeneous  distributions  of  liquid  or  glaciated  precipitation.  The  approach 
which  we  implemented  uses  a  numerical,  Gaussian  quadrature  approach  to  the 
radiative  transfer  equation  (Jin  and  Isaacs,  1985).  The  paper  by  Jin  and 
Isaacs  (1987)  illustrates  applications  of  the  polarized  brightness  temperature 
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calculations  to  the  simulation  of  precipitation  dependent  (i.e.  parameterized 
by  the  rain  rate)  channel  brightness  temperatures  for  the  SSM/I  microwave 
imager . 

3 . 5  Retrieval  Codes 

The  RADTRAN  data  simulation  capability  would  be  greatly  enhanced  by  the 
availability  of  retrieval  codes  to  perform  retrieval  simulations  for  weather 
parameters  of  interest  such  as  temperature  and  moisture  profiles,  cloud  and 
precipitation.  We  have  used  statistical  (Isaacs  and  Deblonde,  1985,  1987)  and 
physical  (Isaacs,  1989)  retrieval  codes  for  this  purpose.  These  have  been 
reviewed  by  Isaacs  (1988b). 

4.  INCORPORATION  OF  ALGORITHM  ENHANCEMENTS 

The  RADTRAN  code  has  been  upgraded  and  modularized  by  incorporating  the 
capabilities  to:  (a)  evaluate  frequency  dependent,  polarized  surface  emis- 
sivity  -  a  menu  driven,  user  selected  surface  type  solution  has  been  imple¬ 
mented  based  on  the  surface  emissivity  submodels  described  in  Section  3.2. 

The  surface  emissivity  supports  evaluation  of  surface  emitted  brightness 
temperatures;  (b)  calculate  scattering  properties  of  precipitation  -  The  sub¬ 
program  for  calculating  precipitation  optical  properties  described  in 
Section  3.3  has  been  incorporated.  This  provides  the  user  with  frequency 
dependent  values  of  extinction  coefficient,  single  scattering  albedo, 
asymmetry  factor,  and  the  angular  scattering  function  for  liquid  and  glaciated 
precipitation.  These  quantities  are  essential  to  performing  the  scalar 
multiple  scattering  calculation  of  brightness  temperatures  in  the  presence  of 
precipitation;  (c)  perform  multiple  scattering  brightness  temperature 
calculations  in  precipitating  conditions:  an  exact  multiple  scattering 
approach  for  fully  polarized  brightness  temperature  calculations  has  been 
included  as  a  user  selectible  option.  This  option  is  applicable  to  the 
simulation  of  data  from  sensors  with  polarization  discrimination;  and 
(d)  perform  statistical  retrievals.  The  specification  of  atmospheric  profiles 
has  also  been  generalized  to  assimilate  user  specified  input  profiles. 

Standard  meteorological  profile  models  (e.g.  U.S.  standard,  tropical,  etc.) 
has  been  augmented  by  the  capability  to  accept  field  data,  radiosonde  and 
other  upper  air  data,  and  user  specified  format  profile  data  as  input. 
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4.1 


The  basic  RADTRAN  algorithm  was  modularized  and  reconfigured  for 
efficiency.  A  main  program  called  RADTRAN  functions  as  the  executive  program 
and  calls  all  other  subroutines.  (Flowcharts  and  complete  user  instructions 
are  given  in  Section  5,  see  Figures  2,  3  and  4.)  RADTRAN  includes  a 
description  of  file  assignments  and  user  input  instructions.  Input  formats 
were  modified  for  consistency  with  the  general  approach  adopted  in  the 
LOWTRAN/FASCODE  environment.  A  major  focus  of  the  algorithm  upgrade  was  the 
integration  of  new  capabilities  to  internally  calculate  surface  emissivities 
from  among  a  variety  of  surface  models,  to  evaluate  both  scattering  and 
attenuation  properties  of  precipitation,  and  to  accomplish  accurate  multiple 
scattering  calculations  in  the  presence  of  precipitation.  Therefore,  user 
input  options  are  provided  to  flag  these  capabilities  on  or  off  and  to  default 
to  the  original  RADTRAN  capability,  i.e.  user  input  of  surface  emissivities 
and  no  multiple  scattering. 

RADTRAN  calls  one  of  five  main  subroutines.  ATMPRF  is  called  to  set  up 
and  print  out  the  user  specified  atmosphere  including  temperature,  water 
vapor,  cloud,  and  rain  profiles.  This  subroutine  controls  the  options  for 
user  supplied,  arbitrary  format  input  model  atmospheres.  Subroutine  ATTEN 
calculates  attenuation  for  the  desired  path  and  also  evaluates  weighting 
functions  if  desired.  RDTRAN  performs  the  radiative  transfer  calculation  in 
non-scattering  atmospheres.  It  incorporates  the  attenuation  calculations  from 
ATTEN  and  the  surface  reflection/emission  either  input  by  the  user  or 
calculated  from  the  surface  emissivity  models.  The  remaining  subroutines 
incorporate  one  of  the  new  capabilities  cited  above  and  interface  with  the 
radiative  transfer  calculation  through  RADTRAN,  depending  on  the  nature  of  the 
options  which  are  toggled  by  the  user. 

4 . 2  Surface  Emissivltv  Enhancement 

Section  3.2  provides  some  background  on  the  internal  surface  model 
capabilities  which  form  the  basis  of  this  enhancement.  The  enhanced  RADTRAN 
implementation  provides  the  user  with  four  options  related  to  the  treatment  of 
the  surface.  These  are  selected  by  toggling  the  flag  ISRFC  in  the  user  input 
file.  The  options  are: 

I 
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ISRFC  -  0 

-  1 


-  2 

-  3 


NO  SURFACE  USED  IN  CALCULATION 
INTERNAL  MODEL  SUPPLIES  SURFACE  TEMPERATURE 
AND  EMISSIVITY 

(NOTE:  FOR  MODELS  1-4,  FREQ  MUST  BE  LESS  THAN  35  GHZ) 
USER  SUPPLIED  SURFACE  TEMPERATURE,  INTERNAL 
MODEL  PROVIDES  EMISSIVITY 
USER  SUPPLIED  SURFACE  TEMPERATURE  AND 
EMISSIVITY 


Option  1  (ISRFC-0)  is  self  explanatory,  i.e.  no  surface  is  used  in  the 
calculation  of  brightness  temperature.  This  implies  no  surface  emission  and 
no  reflection  of  downward  brightness  temperature  contributions.  Option  2 
(ISRFC-1)  selects  one  of  the  model  surfaces  described  briefly  in  Section  3.2 
(see  Table  2) .  The  surface  modeling  subpackage  is  contained  in  subroutine 
SRFMOD  and  its  associated  subroutines.  The  surface  model  selected  is  toggled 
by  the  flag.  Depending  on  the  model  selected  a  variety  of  input  parameters 
may  be  required.  These  are  described  in  the  input  file  description  contained 
in  Section  5.2.  Note  that  subroutine  SRFMOD  can  be  removed  from  the  RADTRAN 
code  and  run  offline  if  it  is  desired  only  to  calculate  polarized  surface 
emissivities .  In  Option  3  (ISRFC-2),  the  user  can  take  advantage  of  the 
ability  of  the  internal  surface  model  subpackage  to  provide  results  for 
varying  surface  related  input  parameters.  The  most  common  surface  related 
parameter  is  the  surface  temperature.  In  Option  2,  the  surface  temperature 
used  in  the  surface  emissivity  calculation  is  chosen  to  match  that  of  the 
lowest  model  layer  of  the  user  selected  input  atmosphere.  In  Option  3,  the 
user  can  select  the  surface  temperature  independently  of  the  atmospheric 
input.  Option  4  (ISRFC-3)  provided  continuity  with  the  basic  RADTRAN  code  and 
requires  that  the  user  supply  both  the  surface  temperature  and  the  surface 
emissivity. 

The  surface  dependent  model  parameters  for  the  various  surface  models 
include  sea  ice  thicknesses  and  scattering  inclusion  size,  snow  depth  and 
wetness,  soil  moisture,  and  vegetative  moisture  and  thickness.  Individual 
users  may  want  to  vary  these  quantities  as  required  for  their  own  simulations 
requirements . 
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4.3  Scattering  Properties  of  Precipitation  Enhancement 

The  basic  RADTRAN  provided  the  capability  to  evaluate  the  attenuation 
properties  of  precipitation,  however,  since  scattering  was  not  included  in  the 
brightness  temperature  calculation,  no  provision  was  made  to  treat  precipita¬ 
tion  scattering  properties.  In  the  enhanced  RADTRAN,  we  have  included  the 
option  of  printing  out  a  table  of  precipitation  scattering  properties  based  on 
the  approximation  described  by  Isaacs  et  al.  (1988).  The  coding  for  this 
function  is  contained  in  subroutine  SCAT.  These  properties  include  the 
absorption  coefficient,  the  extinction  coefficient,  and  the  first  eight 
Legendre  polynomial  coefficients  of  the  expansion  of  the  angular  scattering 
distribution  as  functions  of  frequency,  rainrate,  temperature  phase,  and  size 
distribution  (either  Marshall -Palmer  or  Best).  These  properties  can  be  used 
directly  to  drive  a  simple  scalar  (i.e.  unpolarized)  multiple  scattering 
calculation  of  brightness  temperatures  in  the  presence  of  precipitation  using, 
for  example,  a  discrete  ordinate  multiple  scattering  algorithm  such  as  that 
described  by  Liou  et  al .  (1980).  They  are  also  of  interest  themselves  as  a 
means  to  characterize  potential  scattering  effects  and  contributions  to  the 
brightness  temperature  field  before  multiple  scattering  calculations  are 
performed,  e.g.  by  examining  the  magnitude  of  the  single  scattering  albedo. 
Note  that  these  properties  are  not  used  directly  in  our  multiple  scattering 
enhancement  (Section  4.4)  since  we  require  a  fully  polarized  brightness 
temperature  calculation. 


4 • 4  Polarized.  Multiple  Scattering  Enhancement 

Jin  and  Isaacs  (1985,1987)  describe  a  discrete  ordinate  radiative 
transfer  code  for  the  evaluation  of  polarized  microwave  brightness  tempera¬ 
tures  in  the  presence  of  precipitation.  The  approach  is  applicable  to 
realistic  precipitation  distributions  with  varying  vertical  properties  (i.e., 
inhomogeneity)  such  as  height  dependent  rain  rates  and  rain  phase  (ice/water) 
such  as  those  which  comprise  the  RADTRAN  rain  models.  The  previous  version  of 
RADTRAN  had  no  capability  to  calculate  brightness  temperature  fields  in  the 
presence  of  precipitation  due  to  the  lack  of  a  treatment  for  multiple 
scattering.  (Resulting  brightness  temperatures  including  rain  attenuation 
simply  treated  rain  as  an  equivalent  absorber  thus  failing  to  simulate  several 
important  scattering  mechanisms.) 
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The  enhanced  RADTRAN  incorporates  multiple  scattering  in  a  self 
consistent  manner  and  integrates  it  with  both  the  gas  attenuation  modules  of 
the  original  RADTRAN  (in  subroutine  ATTEN)  and  the  new  surface  modeling 
capabilities  (in  subroutine  SRFMOD) .  The  multiple  scattering  capability  is 
toggled  on  by  the  user  by  the  IMS  input  flag.  Gas  absorption  coefficients 
from  ATTEN  are  input  to  the  multiple  scattering  code  and  combined  with  an 
internal  calculation  of  the  rain  dependent  scattering  and  absorption 
coefficients  to  produce  a  layer  dependent  single  scattering  albedo  profile. 

The  phase  matrix  elements  as  a  function  of  rain  rate  are  also  calculated 
internally.  All  rain  extinction  and  scattering  properties  are  driven  by  the 
user  specified  precipitation  rain  rate/phase  vertical  distribution 
properties.  The  multiple  scattering  code  generates  an  internal  layer  for  each 
layer  of  precipitation  in  the  input  atmosphere.  The  original  Jin  and  Isaacs 
(1987)  algorithm  which  treated  a  two  layer  rain  model  (i.e.  ice  layer  over 
rain  layer  or  two  rain  layers  with  differing  rain  rates)  was  extensively 
modified  to  provide  a  generalized  "n"  layer  calculation.  No  attempt  was  made 
co  reduce  this  number  of  layers  through  layer  merging  procedures.  This  could 
potentially  reduce  computer  time  for  the  multiple  scattering  option. 

The  multiple  scattering  option  runs  interactively  with  the  surface 
modeling  capability.  The  appropriate  surface  emissivity  as  a  function  of 
angle  is  utilized  to  provide  lower  boundary  conditions  to  the  multiple 
scattering  calculation  which  takes  into  account  the  surface  emission,  upward 
atmospheric  emission,  and  reflection  of  downward  atmospheric  emission  incident 
at  the  lower  boundary  of  the  precipitation.  This  boundary  condition  is 
satisfied  for  each  Gauss  point  of  the  discrete  ordinate  calculation. 

4 . 5  Retrieval  Algorithm 

We  have  included  a  generalized  statistical  parameter  inversion  method  for 
use  with  RADTRAN.  The  statistical  retrieval  approach  is  described  in  Isaacs 
and  Deblonde  (1987),  Isaacs  et  al.  (1988),  and  Isaacs  (1987,1988)  applied  to 
the  remote  sensing  of  surface  temperature  and  emissivity,  temperature  profile, 
and  water  vapor  profile.  Jin  and  Isaacs  (1987)  also  applied  the  approach  to 
the  inference  of  rain  rate  using  simulated  scattered  brightness  tempera¬ 
tures.  Isaacs  (1989)  uses  the  statistical  retrieval  as  a  first  guess  for  the 
unified  retrieval  approach  for  DMSP  meteorological  sensors. 
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The  i  ♦-rieval  algorithm  is  not  integrated  with  the  RADTRAN  code,  rather 
it  functions  as  a  separate  entity.  In  application,  the  relationship  between 
the  statistical  retrieval  algorithm  and  RADTRAN  is  a  simple  one.  RADTRAN  is 
used  as  the  forward  problem  calculation  to  generate  synthetic  data  from  an 
ensemble  of  model  atmospheres  or  a  radiosonde  data  set  for  a  desired  sensor 
system.  The  statistical  retrieval  algorithm  reads  in  these  simulated  data  to 
calculate  the  required  statistical  regression  coefficients  for  the  retrieval 
from  one  set  of  simulated  data  and  to  test  the  accuracy  of  the  retrieval  using 
the  regression  statistics  from  a  second  independent  set  of  simulated  data. 

The  procedure  is  described  in  the  above  cited  references. 

5.  DOCUMENTATION  AND  USERS  GUIDE 

5 . 1  Description  of  Enhanced  RADTRAN 

The  original  RADTRAN  algorithm  consisted  of  four  integrated  modules: 

a)  atmospheric  models  (including  rain) 

b)  attenuation  calculation 

c)  calculation  of  weighting  functions 

d)  brightness  temperature  calculation 

All  four  of  the  above  functions  were  performed  in  the  main  body  of  the 
original  RADTRAN  (formerly  named  SKYTMP) .  In  order  to  improve  the  readability 
of  the  code  three  new  subroutines  were  created  to  isolate  computational 
elements  from  the  main  program  which  was  converted  into  an  executive  calling 
routine.  The  new  subroutines  are  (see  Figure  2): 

a)  ATMPRF  which  sets  up  and  prints  out  the  desired  atmospheric 
profile.  ATMPRF  also  alternatively  calls  INPROF  to  read  in  a  user- 
supplied  profile. 

b)  ATTEN  which  performs  the  attenuation  calculation  and  also  calculates 
weighting  functions  if  desired. 

c)  RDTRAN  which  calculates  the  brightness  temperatures  including  the 
contribution  from  the  surface  if  desired. 


16 


The  rotational  lines  for  water  are  now  in  data  statements  and  don't  require 
inputting  unless  the  user  desires  a  new  set. 

The  input  parameters  themselves  were  slightly  modified,  as  can  be  seen 
from  the  user  instructions  (see  Section  5.2).  The  major  changes  however  are 
in  the  three  new  features  that  have  been  added.  These  are:  a)  precipitation 
models  which  provide  the  absorption  and  extinction  coefficients  for  the 
desired  rain  rate,  frequency,  and  temperature;  b)  surface  models  which  provide 
the  surface  emissivity  for  a  variety  of  surfaces  including  ocean,  snow,  sea 
ice,  soil,  and  vegetation;  and  c)  a  multiple  scattering  routine  has  been  added 
so  that  the  scattering  effects  of  precipitating  layers  in  the  atmosphere  can 
be  modeled. 

The  implementation  of  the  precipitation  model  was  accomplished  by  the 
addition  of  three  (3)  new  routines:  a)  SCAT,  b)  BS ,  and  c)  TAB.  SCAT  is  the 
driver  for  the  precipitation  model  and  uses  BS ,  a  binary  search  routine,  and 
absorption  and  extinction  coefficients  for  the  desired  inputs.  The  model 
outputs  are  printed  to  the  main  output  file  RADOUT,  but  are  not  used  directly 
by  RADTRAN . 

The  surface  models  (see  Figure  3)  on  the  other  hand  can  be  used  to 
directly  feed  the  required  emissivities  for  the  desired  look  angles  into  the 
brightness  temperature  calculation.  The  main  driver  SRFMOD  is  used  to  read  in 
the  necessary  inputs,  calls  DWATER  and  DSOIL  to  provide  the  necessary 
dielectric  constants,  calls  EMIS  to  produce  the  emissivities  at  the  gauss 
points,  and  then  calls  AITINT  to  interpolate  the  emissivities  to  the  desired 
angles.  EMIS  is  the  primary  routine  in  which  the  code  decides  based  on  user 
input  how  to  calculate  the  emissivities.  For  calm  ocean  the  program  utilizes 
a  model  of  a  dielectric  slab,  for  rough  ocean,  sea  ice,  or  snow,  a  model  of 
random  discrete  scatterers  is  used,  and  for  vegetation  or  soil,  a  continuous 
random  medium  model  is  used. 

The  multiple  scattering  routine  (see  Figure  4)  provides  another  alter¬ 
native  to  the  standard  RADTRAN  run.  It  allows  the  user  to  fully  model  the 
effects  that  precipitating  layers  have  on  the  observed  radiance.  The  driver 
routine  MLTSCT  sets  up  the  layers  containing  precipitation  for  the  model,  and 
calculates  the  boundary  conditions  using  the  routine  RDTNMS .  The  polarized 
multiple  scattered  radiance  is  then  calculated  in  MSPTRT.  This  radiance  which 
is  calculated  at  the  six  (6)  Gauss  points  is  then  interpolated  by  AITINT. 
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SRFMOD 


Figure  3.  Surface  model  subroutines  flowchart. 


MLTSCT 
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Figure  4.  Multiple  scattering  model  subroutines  flowchart 


The  program  structure  and  subroutine  hierarchy  of  enhanced  RADTRAN  is 
outlined  in  Table  3,  which  also  provides  brief  descriptions  of  all  subroutines 
and  functions.  What  is  not  obvious  however,  is  the  scope  of  the  changes  that 
have  been  made.  RADTRAN  as  it  previously  existed  has  approximately  1500  lines 
of  code.  The  new  enhanced  RADTRAN  is  approximately  7500  lines  of  code. 
Unfortunately  along  with  the  increase  in  size  of  the  code  came  an  increase  in 
the  storage  required.  Currently  the  code  will  run  unsegmented  on  AFGL's  CYBER 
180-860,  it  does  however  use  approximately  300K  of  extended  core. 

Table  3.  Subroutine  Hierarchical  Structure  for  the  Enhanced  RADTRAN 


PROGRAM  RADTRAN 
|  SUBROUTINE  SCAT 
|  |  SUBROUTINE  BS 
|  |  FUNCTION  TAB 

I  I  I 

|  SUBROUTINE  ATMPRF 
|  |  SUBROUTINE  INPROF 
|  |  FUNCTION  RHOSAT 
|  SUBROUTINE  SRFMOD 
|  |  SUBROUTINE  DWATER 
|  j  SUBROUTINE  DSOIL 
|  |  |  SUBROUTINE  DIE 
|  |  SUBROUTINE  EMIS 
|  |  |  SUBROUTINE  DIE 
|  |  |  SUBROUTINE  DRYDS 

I  I  I  I 

|  |  |  SUBROUTINE  WETDS 

II  I  I 

|  |  |  |  SUBROUTINE  ZANLYT 

I  I  I  I  I 

|  |  |  |  |  SUBROUTINE  UERTST 
|  |  |  |  ||  SUBROUTINE  USPKD 

|  |  |  |  |  |  SUBROUTINE  UGETIO 

1111  FUNCTION  CFUN 

I  I  I  II 

|  |  |  SUBROUTINE  RANDOM 
I  I  I  I 

|  |  |  |  SUBROUTINE  VEGE 

I  I  I  II 

|  |  |  |  SUBROUTINE  SOIL 

I  I  I  II 

|  |  |  |  SUBROUTINE  REF 
1111  SUBROUTINE  EBIS 
|  |  |  |  |  SUBROUTINE  QSF 


Main  Program 

Driver  for  Precipitation  Model 
Binary  Search  Routine 
Table  of  Absorption  and 
Extinction  Coefficients 
Sets  up  the  Atmospheric  Profile 
Reads  in  User  Supplied  Profile 
Computes  Water  Vapor  Density 
Driver  for  Surface  Models 
Dielectric  Constant  for  Water 
Dielectric  Constant  for  Soil 
Water  and  Ice  Dielectric  Permittivities 
Driver  for  Emissivity  Calculation 
Water  and  Ice  Dielectric  Permittivities 
Calculate  Emissivities  for  a  Dry 

Random  Discrete  Scattering  Medium 
Calculate  Emissivities  for  a  Wet 

Random  Discrete  Scattering  Medium 
IMSL  -  Zeroes  of  an  Analytic 
Complex  Function 
IMSL  -  Prints  Error  Condition 
IMSL  -  Required  by  UERTST 
IMSL  -  Obtains  I/O  Units 
Algebraic  Function 
X3+C1*X2+C2*X+C3 
Calculate  Emissivities  for 
Continuous  Random  Medium 
Mean  Dielectric  Constant  and 
Variance  for  Vegetation 
Mean  Dielectric  Constant  and 
Variance  for  Soil 
Calculates  the  Reflectivities 
Compute  Scattering  Phase  Functions 
Compute  Vector  of  Integral  Values 
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Table  3.  Subroutine  Hierarchical  Structure  for  the  Enhanced  RADTRAN 
(continued) 


|  |  |  SUBROUTINE  EIGRF 

Ml! 

Mil  SUBROUTINE  EBALAF 
Mil  SUBROUTINE  EHESSF 
|  |  |  |  SUBROUTINE  EHBCKF 
|  |  |  |  SUBROUTINE  EQRH3F 
jjlll  SUBROUTINE  UERTST 
|  |  |  |  |  |  SUBROUTINE  USPKD 
||  I  II  I  SUBROUTINE  UGETIO 
|  |  |  I  SUBROUTINE  EBBCKF 
|  |  |  I  SUBROUTINE  UERTST 
jjlll  SUBROUTINE  USPKD 
|  |  |  |  |  SUBROUTINE  UGETIO 
|  |  |  SUBROUTINE  SIMQ 
j  j  SUBROUTINE  REF 
|  SUBROUTINE  AITINT 
SUBROUTINE  RDH20 
SUBROUTINE  PLTID3 
SUBROUTINE  PLOT 
SUBROUTINE  ATTEN 
|  FUNCTION  GAM02 
|  FUNCTION  GAMH20 
|  |  SUBROUTINE  ALPHAG 
|  |  SUBROUTINE  ALPHAS 
|  FUNCTION  GAMFOG 

I  I 

|  |  SUBROUTINE  DPH20 
|  |  SUBROUTINE  INDEX 
j  j  |  SUBROUTINE  DEBYE 
|  |  |  FUNCTION  DOP 
|  |  |  FUNCTION  AB 
|  FUNCTION  GMRAIN 

I  I 

|  |  FUNCTION  AITK 
SUBROUTINE  MLTSCT 
|  SUBROUTINE  RDTNMS 

II 

|  SUBROUTINE  MSPTRT 
|  |  SUBROUTINE  TSF 
|  |  SUBROUTINE  DWATER 
j  |  SUBROUTINE  DIE 
j  j  SUBROUTINE  QSF 
|  |  SUBROUTINE  EIGRF 

I  I  I 

|  |  |  SUBROUTINE  EBALAF 

|  |  j  SUBROUTINE  EHESSF 

j  j  j  SUBROUTINE  EHBCKF 

|  |  |  SUBROUTINE  EQRH3F 


IMSL  -  Eigenvalues  and  Eigenvectors 
of  a  Real  General  Matrix 
IMSL  -  Required  by  EIGRF 
IMSL  -  Required  by  EIGRF 
IMSL  -  Required  by  EIGRF 
IMSL  -  Required  by  EIGRF 
IMSL  -  Prints  Error  Condition 
IMSL  -  Required  by  UERTST 
IMSL  -  Obtains  I/O  Units 
IMSL  -  Required  by  EIGRF 
IMSL  -  Prints  Error  Condition 
IMSL  -  Required  by  UERTST 
IMSL  -  Obtains  I/O  Units 
Solve  Matrix  AX-B 
Calculates  the  Reflectivities 
Aitken's  Iterated  Interpolation 
Reads  in  H20  Rotational  Lines 
AFGL  Plot  Header  Routine 
Various  Plot  Functions 
Compute  Atmospheric  Attenuation 
Compute  Oxygen  Attenuation 
Driver  for  Computing  Water  Attenuation 
Computes  Water  Attenuation  (>300  GHZ) 
Computes  Water  Attenuation  (5<f<300  GHZ) 
Computes  Attenuation  for  Condensed  Water 
in  Cloud  or  Fog 

Computes  Dielectric  Properties  of  Water 

Computes  Refractive  Index 

Required  by  INDEX 

Required  by  INDEX 

Required  by  INDEX 

Computes  Attenuation  for  Condensed  Water 
in  form  of  Rain 
Required  by  GMRAIN 

Driver  for  Multiple  Scattering  Routine 
Computes  Boundary  Conditions  and 

Contributions  from  non-scattering  layer 
Multiple  Scattering  Routine 
Computes  TNM  and  SNM  Functions 
Dielectric  Constant  for  Water 
Water  and  Ice  Dielectric  Permittivities 
Compute  Vector  of  Integral  Values 
IMSL  -  Eigenvalues  and  Eigenvectors 
of  a  Real  General  Matrix 
IMSL  -  Required  by  EIGRF 
IMSL  -  Required  by  EIGRF 
IMSL  -  Required  by  EIGRF 
IMSL  -  Required  by  EIGRF 
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Table  3.  Subroutine  Hierarchical  Structure  for  the  Enhanced  RADTRAN 
(continued) 


|  |  |  |  SUBROUTINE  UERTST 
jjjjl  SUBROUTINE  USPKD 
Mill  SUBROUTINE  UGETIO 
j  |  |  SUBROUTINE  EBBCKF 
|  j  j  SUBROUTINE  UERTST 
||  I  I  SUBROUTINE  USPKD 
||  I  I  SUBROUTINE  UGETIO 
|  |  SUBROUTINE  SIMQ 
j  SUBROUTINE  AITINT 
SUBROUTINE  RDTRAN 

I 

SUBROUTINE  DRAW 
|  SUBROUTINE  AXIS 
|  SUBROUTINE  NUMBER 
|  SUBROUTINE  LINE 
|  SUBROUTINE  SYMB 
j  |  SUBROUTINE  SYMBOL 
j  SUBROUTINE  WHERE 
j  SUBROUTINE  PLOT 
|  SUBROUTINE  FRAME 
j  SUBROUTINE  ENDPLT 


IMSL  -  Prints  Error  Condition 
IMSL  -  Required  by  UERTST 
IMSL  -  Obtains  I/O  Units 
IMSL  -  Required  by  EIGRF 
IMSL  -  Prints  Error  Condition 
IMSL  -  Required  by  UERTST 
IMSL  -  Obtains  I/O  Units 
Solve  Matrix  AX-B 
Aitken's  Iterated  Interpolation 
Computes  Brightness  Temperature 

for  the  non-Multiple  Scattering  Cases 
Driver  for  Plotting  Routines 
Draws  Axis 
Draws  Number 
Draws  a  Line 
Sets  up  call  to  SYMBOL 
Draws  a  Symbol 
Returns  Pen  Location 
Various  Plot  Functions 
Draws  Frame 

AFGL  Plot  Terminating  Routine 


5 . 2  Enhanced  RADTRAN  User  Instructions 

Although  the  enhanced  RADTRAN  code  now  contained  some  sophisticated 
options  every  attempt  was  made  to  keep  the  user  interlace  as  simple  as 
possible.  This  meant  continuity  with  the  old  RADTRAN  input  format  and  the 
introduction  of  LOWTRAN  type  input  and  file  assignments.  Table  4  lists  the 
file  assignments  for  RADTRAN.  Table  5  summarizes  the  input  instructions  for 
enhanced  RADTRAN  as  they  are  included  in  the  code  itself. 

Table  4.  File  Assignments  for  RADTRAN 


TAPES  RADTRAN  INPUT  FILE  (SKYPAR) 

TAPE6  RADTRAN  OUTPUT  FILE  (RADOUT) 

TAPE7  TABULAR  RADIANCE  OUTPUT  FILE  (RADTAB) 

TAPE8  SURFACE  MODEL  OUTPUT  FILE  (EMIDAT) 

TAPE9  BINARY  OUTPUT  FILE  -  USED  TO  STORE 

WEIGHTING  FUNCTIONS  (RADPLT) 
TAPE10  SPECTRAL  RANGE  OF  FREQUENCIES  (RADSPC) 

TAPE20  WATER  VAPOR  SPECTRAL  LWE  PARAMETER  INPUT  FILE  (H20R0T) 

TAPE39  AFGL  PLOT  FILE 
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Table  5.  Input  Instructions  for  RADTRAN 


CARD  1 :  TITLE 
FORMAT  (A80) 

TITLE  IS  USED  AS  A  HEADER  RECORD  FOR  OUTPUT  PURPOSES 


CARD  2:  IRDTRN , IATM, ISRFC , IMS , I PLOT , NANGLE , MODE , ITABLE, 

IATTEN , IATINC , IWF, IWFINC, IRDH20 , NOPR , TGRND 
FORMAT  (14(3X, 12) , F10 . 5) 

IRDTRN-  0  PROFILE  AND  SURFACE  MODELS  SET  UP  (RADTRAN  NOT  RUN) 

-  1  FOR  RADTRAN  RUN 

-  2  FOR  RUN  OF  PRECIPITATION  MODEL  (RADTRAN  NOT  RUN) 

IATM  -  0  FOR  USER  SUPPLIED  ATMOSPHERE 

-  1  FOR  STANDARD  MODEL 

ISRFC  -  0  NO  SURFACE  USED  IN  CALCULATION 

-  1  INTERNAL  MODEL  SUPPLIES  SURFACE  TEMPERATURE 

AND  EMISSIVITY 

(NOTE:  FOR  MODELS  1-4,  FREQ  MUST  BE  LESS  THAN  35  GHZ) 

-  2  USER  SUPPLIED  SURFACE  TEMPERATURE,  INTERNAL 

MODEL  PROVIDES  EMISSIVITY 

-  3  USER  SUPPLIED  SURFACE  TEMPERATURE  AND 

EMISSIVITY 

IMS  -  0  NO  MULTIPLE  SCATTERING 

-  1  MULTIPLE  SCATTERING  RUN 

I PLOT  -  0  NO  PLOTTING 

-  1  PLOT 

NANGLE  IS  THE  NUMBER  OF  ANGLES  TO  BE  INPUTED  (MAX-10) 

MODE  DETERMINES  HOW  FREQUENCIES  ARE  INPUTED 

-  0  READ  IN  BEGINNING  AND  ENDING  FREQUENCY  AND  INCREMENT 

-  1  READ  IN  INDIVIDUAL  FREQUENCIES 

ITABLE  DETERMINES  IF  TABULAR  OUTPUT  IS  WRITTED  TO  RADTAB 

-  0  FOR  NO  TABLE  PRINTED 

-  1  FOR  TABLE  PRINTED  TO  RADTAB 

IATTEN  DETERMINES  IF  ATTENUATIONS  ARE  PRINTED 

-  0  FOR  NO  ATTENUATION  PRINTED 

-  1  FOR  ATTENUATION  PRINTED 

IATINC  IS  THE  INCREMENT  FOR  WHICH  ATTENUATIONS  ARE  PRINTED 

-  1  EVERY  ATTENUATION  PRINTED 

-  N  FOR  EVERY  NTH  ATTENUATION  PRINTED 
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Table  5.  Input  Instructions  for  RADTRAN  (continued) 


IWF  DETERMINES  IF  WEIGHTING  FUNCTIONS  ARE  PRINTED 

-  0  FOR  NO  WEIGHTING  FUNCTIONS  PRINTED 

-  1  FOR  WEIGHTING  FUNCTIONS  PRINTED 

IWFINC  IS  THE  INCREMENT  FOR  WHICH  WEIGHTING  FUNCTIONS  ARE  PRINTED 

-  1  EVERY  WEIGHTING  FUNCTION  PRINTED 

-  N  FOR  EVERY  NTH  WEIGHTING  FUNCTION  PRINTED 

IRDH20  IS  THE  FLAG  FOR  USER  INPUTED  WATER-VAPOR  ROTATION  SPECTRAL 
LINE  PARAMETERS 

-  0  USES  INTERNAL  LINES 

-  1  READS  IN  LINES  FROM  FILE  H20R0T 

NOPR  IS  THE  FLAG  TO  RESTRICT  THE  SIZE  OF  THE  PRINTOUT 

-  0  FULL  PRINTOUT 

-  1  LIMITED  PRINTOUT 

TGRND  IS  THE  SURFACE  TEMPERATURE  (ISRFC=2,3) 

(SET  TO  T(l)  FOR  TGRND-O.O) 


CARD  3:  DIST , PHASE, FREQ, RNRATE, TEMP  (IRDTRN-2) 

FORMAT  (3F10. 3 , 215) 

DIST  DETERMINES  THE  DROP  SIZE  DISTRIBUTION 

-  1  FOR  MARSHALL- PALMER 

-  2  FOR  BEST 

-  -1  WILL  TERMINATE  PRECIPITATION  MODEL  RUN 

PHASE  DETERMINES  THE  PHASE  OF  THE  PRECIPITATION 
=  1  FOR  WATER 
-  2  FOR  ICE 

FREQ  IS  THE  FREQUENCY  IN  GHZ  (19.35  -  231.0  GHZ) 

RNRATE  IS  THE  RAINFALL  RATE  IN  MM/HR  (0.0  -  50.0  MM/HR) 
TEMP  IS  THE  TEMPERATURE  IN  DEGREES  K  (263.0  -  283.0  K) 
****  REPEAT  CARD  3  **** 


CARD  4.1.1:  ICNTRL , M , ATITLE  (I ATM-0) 
FORMAT  (2(3X,I2) ,A70) 
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Table  5.  Input  Instructions  for  RADTRAN  (continued) 


ICNTRL  DETERMINES  HOW  PROFILE  IS  TO  BE  INPUTED 

-  0  READ  IN  PROFILE  IN  TABULAR  FORM 

-  1  READ  IN  PROFILE  IN  COMPACT  FORM 

-  2  READ  IN  PROFILE  IN  COMPACT  FORM  (NO  RAIN) 

-  3  READ  IN  PROFILE  IN  COMPACT  FORM  (NO  CLOUDS,  NO  RAIN) 

M  IS  THE  NUMBER  OF  LEVELS  TO  BE  READ  IN 

ATITLE  IS  A  HEADER  TO  IDENTIFY  THE  PROFILE  INPUTED 


CARD  4.1.2:  FMT  (ICNTRL-0,1) 

FORMAT  (A80) 

FMT  IS  THE  FORMAT  TO  BE  USED  FOR  INPUTING  THE  PROFILE 


CARD  4.1.3:  H.P.T.RH.CLOUD.RAIN  (ICNTRL-0) 
FORMAT  (INPUT  FROM  CARD  4.1.2) 

H  IS  THE  HEIGHT  IN  KM 

P  IS  THE  PRESSURE  IN  MB 

T  IS  THE  HEIGHT  IN  KM 

RH  IS  THE  RELATIVE  HUMIDITY 

CLOUD  IS  THE  CLOUD  AMOUNT 

RAIN  IS  THE  RAIN  RATE  (MM/HR) 


CARD  4.2:  MOD , MHUMID , MCLOUD , MRAIN  (IATM-1) 
FORMAT  (4(3X,I2)) 


MOD 


-  1  TROPICAL  MODEL 

-  2  MID-LATITUDE  SUMMER  MODEL 

-  3  MID-LATITUDE  WINTER  MODEL 

-  4  SUB-ARCTIC  SUMMER  MODEL 

-  5  SUB-ARCTIC  WINTER  MODEL 

-  6  U.S.  STANDARD  (1962)  MODEL 


MHUMID-  1 

-  2 

-  3 

-  4 


TROPICAL  MODEL 
MID -LATITUDE  SUMMER  MODEL 
MID- LATITUDE  WINTER  MODEL 
SUB-ARCTIC  SUMMER  MODEL 
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Table  5.  User  Instructions  for  RADTRAN  (continued) 


-  5  SUB-ARCTIC  WINTER  MODEL 
-6  U.S.  STANDARD  (1962)  MODEL 

-  7  TYPICAL  HUMIDITIES  IN  RAINY  ATMOSPHERE 

-  8  TYPICAL  HUMIDITIES  IN  CLOUD  ONLY  ATMOSPHERE 


MCLOUD-  1 
-  2 

-  3 

-  4 

-  5 

-  6 


NO  CLOUDS 

STRATUS/STRATO  CUMULUS  0.15  G/M3  0. 5-2.0  KM 
CUMULUS  1.0  G/M3  1.0-3. 5  KM 
ALTO  STRATUS  0.4  G/M3  2. 5-3.0  KM 
STRATUS  CUMULUS  0.55  G/M3  0.5-1.5KM 
NIMBOSTRATUS  0.61  G/M3  0.5-1.0KM 


MRAIN  -  1  NO  RAIN 

-  2  DRIZZLE  2 MM/HR  AT  SFC 

-  3  LIGHT  RAIN  5 MM/HR  AT  SFC 

-  4  STEADY  RAIN  12.5MM/HR  AT  SFC 

-  5  SUMMER  CUMULUS  RAIN  15,0  MM/HR  AT  SFC 

-  6  HEAVY  RAIN  20MM/HR  AT  SFC 


CARD  5.1.1:  IGRND, IEMTAB , TEMPI , TEMP2  (ISRFC-1,2) 
FORMAT  (3X,I2,2F10.5) 

IGRND  IS  THE  SURFACE  MODEL  DESIRED 

IEMTAB  IS  THE  FLAG  TO  PRINT  A  TABLE  OF  EMISSIVITES 
TO  THE  FILE  EMI DAT 

-  0  NO  TABLE  PRINTED 

-  1  EMISSIVITY  TABLE  PRINTED 

TEMPI  IS  THE  TEMPERATURE  OF  THE  UPPER  SURFACE  LAYER 
(SET  TO  TGRND  FOR  TEMPI  -  0.0) 

TEMP2  IS  THE  TEMPERATURE  OF  THE  LOWER  SURFACE  LAYER 
(SET  TO  TGRND  FOR  TEMP2  -  0.0) 


CARD  5.1.2:  A , FRAC . EB , E2DRY , DEPDRY  ( IGRND- 1 ,2,3) 
FORMAT  (7F10.6) 

**  ROUGH  OCEAN,  SEA  ICE,  OR  DRY  SNOWPACK  ** 

A  IS  THE  SCATTERER  SIZE 
FRAC  IS  THE  FRACTION  VOLUME 
EB  IS  THE  LAYER  PERMITTIVITY 
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Table  5.  Input  Instructions  for  RADTRAN  (continued) 


E2DRY  IS  THE  PERMITTIVITY  FOR  THE  SURFACE 
DEPDRY  IS  THE  DEPTH  OF  THE  LAYER 


CARD  5.1.3:  Al , FS1 ,A2 , FS2 ,E2WET , DEPWET  (IGRND-4) 
FORMAT  (7F10.6) 

**  WET  SNOWPACK  ** 

Al  IS  THE  SCATTERER  SIZE  FOR  THE  UPPER  LAYER 
FS1  IS  THE  TOTAL  FRACTION  FOR  THE  UPPER  LAYER 
A2  IS  THE  SCATTERER  SIZE  FOR  THE  LOWER  LAYER 
FS2  IS  THE  TOTAL  FRACTION  FOR  THE  LOWER  LAYER 
E2WET  IS  THE  PERMITTIVITY  FOR  LOWER  LAYER 
DEPWET  IS  THE  DEPTH  OF  THE  LAYER 


CARD  5.1.4:  LP , LZ , MV , E2RAN , DEPRAN , H , DMV  ( IGRND-5 ,6,7) 
FORMAT  (8F10.6) 

**  VEGETATION,  WET  SOIL,  AND  RANDOM  MEDIUM  ** 

LP  IS  THE  CORRELATION  LENGTH 
LZ  IS  CORRELATION  OF  Z- EXPONATIOAL 
MV  IS  VOLUMETRIC  MOISTURE  IN  VEGETATION 
E2RAN  IS  THE  DIELECTRIC  CONSTANT 
DEPRAN  IS  THE  DEPTH  OF  THE  LAYER 
H  IS  THE  ROUGHNESS 

DMV  IS  THE  VOLUMETRIC  MOISTURE  USED  IN  CALCULATING  E2RAN 


CARD  5.1.5:  W , VB , VBW , BULK  (IGRND-5) 
FORMAT  (5F10.6) 

**  VEGETATION  ** 

W  IS  THE  VEGETATION  FRACTION  IN  LAYER 
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Table  5.  Input  Instructions  for  RADTRAN  (continued) 


VB  IS  THE  DRY  VEGETATION  IN  VEGETATION 
VBW  IS  THE  BOUNDED  WATER  FRACTION  IN  MV 
BULK  DIELECTRIC  CONSTANT 

CARD  5.1.6:  WT , GA , PO , ROCK  (IGRND-6) 
FORMAT  (5F10.6) 

**  WET  SOIL  ** 

WT  IS  THE  TRANS ION  WETNESS 

GA  IS  THE  MIXING  PARAMETER 

PO  IS  THE  POROCITY 

ROCK  IS  THE  ROCK  PERMITTIVITY 


CARD  5.1.7:  EM, DEL  (IGRND-7) 

FORMAT  (3F10.6) 

**  RANDOM  MEDIUM  ** 

EM  IS  THE  AVERAGE  DIELECTRIC  CONSTANT  OF  THE  RANDOM  MEDIUM 
DEL  IS  THE  VARIANCE 

**  NOTE  :  THESE  ARE  FREQUENCY  DEPENDENT,  SO  FOR  MODE-1 
MUST  ALSO  READ  IN  ON  CARD  9.2 


CARD  5.2:  EMISH(l) , EMISH(2) . EMISH(IO)  (ISRFC-3) 

FORMAT  (10F8.4) 

EMISH(I)  IS  THE  HORIZONTAL  EMISSIVITY  CORRESPONDING 
TO  ANGLE(I) 


CARD  5.3:  EMISV(l) , EMISV(2) . EMISV(IO)  (ISRFC-3) 

FORMAT  (10F8.4) 

EMISV(I)  IS  THE  VERTICAL  EMISSIVITY  CORRESPONDING 
TO  ANGLE(I) 


Table  5.  Input  Instructions  for  RADTRAN  (continued) 

CARD  6.1:  EMHMS  (IMS-1  AND  ISFRC-3) 

FORMAT  (6F10.6) 

EMHMS  IS  THE  HORIZONTAL  EMISSIVITY  CORRESPONDING  TO  EACH  OF 
THE  SIX  GAUSS  POINTS  USED  BY  MULTIPLE  SCATTERING 

GAUSS  ANGLES  ARE  11.02,  25.30,  39.65,  54.03,  68.42,  6.  82.81 


CARD  6.2:  EMVMS  (IMS-1  AND  ISFRC-3) 

FORMAT  (6F10.6) 

EMVMS  IS  THE  VERTICAL  EMISSIVITY  CORRESPONDING  TO  EACH  OF 
THE  SIX  GAUSS  POINTS  USED  BY  MULTIPLE  SCATTERING 

GAUSS  ANGLES  ARE  11.02,  25.30,  39.65,  54.03,  68.42,  &  82.81 


CARD  7:  NAME  (IPLOT-1) 

FORMAT  (A30) 

NAME  IS  THE  CHARACTER  RECORD  USED  TO  IDENTIFY  PLOTS 


CARD  8:  ANGLE(l) ,ANGLE(2) , . . . .ANGLE (10) 
FORMAT  (10F8.4) 

ANGLE(I)  IS  THE  LOOK  ANGLE  (IN  DEGREES) 


CARD  9.1:  V1.V2.DV  (ONLY  FOR  MODE-O) 
FORMAT  (3F10.5) 

VI  IS  THE  BEGINNING  FREQUENCY 

V2  IS  THE  ENDING  FREQUENCY 

DV  IS  THE  FREQUENCY  INCREMENT 


CARD  9.2:  VI ,V2 , V3 , E2 ,DELF  (ONLY  FOR  MODE-1) 
FORMAT  (6F10.5) 

VI  IS  THE  CENTER  FREQUENCY 
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Table  5.  Input  Instructions  for  RADTRAN  (continued) 


V2  IS  THE  FIRST  SIDEBAND 
V3  IS  THE  SECOND  SIDEBAND 

E2  IS  THE  PERMITTIVITY  (ISRFC-1,2  AND  IGRND-1,5,6) 
(PROGRAM  WILL  CALCULATE  FOR  E2-(0. ,0.)) 

E2  IS  (EM)  THE  AVERAGE  DIELECTRIC  CONSTANT  FOR  THE 
RANDOM  MEDIUM  (ISFRC-1,2  AND  IGRND-7) 

DELF  IS  THE  VARIANCE  FOR  EM  (IGRND-7) 


CARD  10:  IRPT 
FORMAT  (3X, 12) 

IRPT  DETERMINES  NEXT  PROGRAM  OPERATION 

-  0  PROGRAM  TERMINATES 

-  1  PROGRAM  RERUN  READING  ALL  INPUT  CARDS 

-  2  PROGRAM  RERUN  READING  IN  NEW  HEADER  AND 

OTHER  INPUTS  STARTING  WITH  ATMOSPHERE 

-  3  PROGRAM  RERUN  READING  IN  NEW  HEADER  AND 

OTHER  INPUTS  STARTING  WITH  SURFACE  TERMS 
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In  order  to  demonstrate  some  of  the  new  features  of  the  RADTRAN  program, 
we  have  generated  four  different  test  cases.  In  test  case  1  we  have  exercised 
the  precipitation  model.  The  precipitation  model  option  (IRDTRN-2)  on  CARD  2, 
allows  the  user  to  obtain  the  absorption  and  extinction  coefficients  along 
with  the  first  eight  legendre  coefficients  as  outlined  in  the  user 
instructions  (see  CARD  3  in  Table  5).  In  the  input  file  for  test  case  1  (see 
Table  6)  we  have  chosen  three  frequencies  (19.35,  37.00,  and  85.50  GHZ),  four 
rainrates  (25.0,  10.0,  5.0,  and  1.0  MM/HR),  and  two  temperatures  (270.0  and 
280.0  K) .  The  output  (see  Table  7)  shows  the  desired  quantities  for  each  of 
the  input  parameters  selected. 

For  the  second  test  case,  the  surface  model  option  is  exercised.  The 
internal  surface  model  option  (ISRFC-1,2)  on  CARD  2,  allows  the  user  to  tailor 
the  surface  parameters  used  by  RADTRAN  to  reflect  several  different  physical 
ground  conditions.  In  the  input  file  for  test  case  2  (see  Table  8),  we  have 
chosen  four  different  surfaces  to  evaluate.  These  are  a)  CALM  OCEAN,  b)  WET 
SNOWPACK,  c)  VEGETATION,  and  d)  WET  SOIL.  Although  we  allowed  most  of  the 
model  parameters  to  use  the  defaults,  a  large  number  of  possible  permutations 
exist  for  the  user  to  chose  from  (see  CARDS  5.1.1  -  5.1.7  in  Table  5).  This 
test  case  produces  three  output  files.  The  main  RADTRAN  output  (see  Table  9) 
is  placed  on  file  RADOUT. 

In  the  first  section  of  this  output  file  is  printed  the  atmospheric 
profile  used.  In  this  case,  a  tropical  model  was  chosen.  Following  the 
atmospheric  profile  is  a  listing  of  the  t*cr *'*i*por  rc tut i wiia  1  spectral  line 
parameters.  These  are  stored  internally,  but  can  be  inputted  by  the  user  if 
desired.  The  next  section  shows  the  attenuations  at  every  level.  This  table 
is  optionally  printed  by  changing  the  IATTEN  flag  on  CARD  2.  Following  the 
attenuations  is  a  table  of  weighting  functions  for  both  ground- to-space  and 
space -to -ground.  The  weighting  functions  can  also  be  optionally  chosen  by 
changing  the  IWF  flag  on  CARD  2.  Following  the  weighting  functions  is  the 
result  for  the  first  frequency  (10.6  GHZ).  This  provides  the  brightness 
temperature  looking  up,  and  the  polarized  brightness  temperature  looking 
down.  Also  provided  is  the  polarized  emissivities  used  by  the  routine.  In 
this  case,  the  emissivities  printed  are  those  obtained  from  the  surface 
model.  The  next  table  contains  the  attenuations  for  the  second  frequency 
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(18.0  GHZ),  followed  by  its  weighting  functions  and  the  brightness  temperature 
results . 

The  next  section  begins  the  RADTRAN  run  for  the  WET  SNOW  case.  The  user 
should  note  that  CARD  10,  allows  the  repeated  runs  of  RADTRAN  from  one  input 
file.  The  rest  of  the  output  file  contains  the  attenuations  and  brightness 
temperature  for  the  VEGETATION  and  WET  SOIL  runs. 

When  running  the  surface  models,  a  more  extensive  output  is  usually 
desired.  This  can  be  obtained  by  setting  IEMTAB-1  on  CARD  5.1.1.  This  will 
produce  the  additional  output  file  EMIDAT  (see  Table  10).  This  output  file 
shows  the  inputs  used  by  the  surface  models  whether  user  selected  or  defaults, 
along  with  the  calculated  emiss ivities .  An  additional  output  which  can  be 
selected  (ITABLE-1  on  CARD  2)  is  printed  on  file  RADTAB  (see  Table  11).  This 
option  provides  a  concise  tabulated  output  file  in  addition  to  the  main 
output . 

For  the  third  test  case,  the  multiple  scattering  option  is  exercised. 

This  option  is  chosen  by  selecting  IMS-1  on  CARD  2.  For  this  test  case  (see 
Table  12)  we  have  chosen  a  moist  vegetation  surface  and  a  tropical  model.  In 
addition  we  have  chosen  two  different  rainrates,  a  light  rain,  and  a  summer 
cumulus  rain.  Both  multiple  scattered  and  no  multiple  scattered  results  are 
desired.  The  only  difference  of  note  in  the  output  files,  is  that  for 
multiple  scattered  cases,  no  upward  brightness  temperatures  are  currently 
calculated.  The  main  output  for  the  test  case  3  run  is  provided  in 
Table  13.  The  surface  model  output  (Table  14)  and  the  tabulated  output 
(Table  15)  are  also  provided.  As  can  be  seen  from  the  output,  there  is  a 
significant  difference  between  the  calculated  brightness  temperatures  from  the 
multiple  scattered  cases  and  the  non  multiple  scattered  cases. 

The  final  test  case  demonstrates  the  MODE-O  option  on  CARD  2.  This 
option  allows  the  user  to  select  a  starting  and  ending  frequency  along  with  an 
increment  to  obtain  the  output  over  a  spectral  range  of  frequencies.  For  this 
case  (see  Table  16)  we  chose  the  interval  182.0-184.0  GHZ,  incrementing  every 
0.1  GHZ.  For  computational  efficiency  the  emissivities  were  set  to  1.0.  The 
main  output  file,  RADOUT  (see  Table  17),  can  be  quite  large  especially  if  the 
attenuations  and  weighting  functions  are  desired.  For  this  case  we  did  not 
turn  on  those  print  options.  In  this  case,  the  advantages  of  the  RADTAB 
output  file  (see  Table  18)  are  readily  apparent.  If  the  user  does  desire  the 


total  attenuations  for  each  frequency,  the  IATTEN-1  flag  produces  an 
additional  output  file,  RADSPC  (see  Table  19).  This  provides  a  concise 
tabulated  file  of  total  attenuations  at  every  frequency.  Note,  however, 

RADOUT  will  be  significantly  larger  with  IATTEN-1. 

An  additional  note  about  cpu  usage.  Each  of  these  test  cases  was  run  on 
AFGL's  CYBER  180-860  and  the  run  time  recorded.  For  case  1,  it  took  .273 
seconds;  for  case  2,  106.524  seconds;  for  case  3,  576.093  second;  and  for  case 
4,  4.835  seconds. 
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Table  6 .  Input  ftfr  Precipitation  Model  Test  Case 


o 

O 

O 

O 

w 

O 

o 

o 

O 

O 

O 

o 

o 

O 

O 

o 

O 

o 

o 

O 

O 

O 

O 

0 

O 

O 

O 

O 

0 

O 

O 

0 

O 

O 

0 

O 

O 

0 

O 

O 

O 

0 

0 

0 

0 

0 

0 

O 

O 

o 

O 

O 

O 

O 

o 

o 

o 

O 

O 

o 

o 

O 

O 

o 

O 

o 

o 

O 

O 

O 

O 

0 

O 

O 

O 

O 

0 

O 

O 

0 

O 

O 

0 

O 

O 

0 

O 

O 

O 

0 

0 

0 

0 

0 

0 

O 

oo 

oo 

CO 

ao 

p-» 

r» 

r-v 

oo 

oo 

oo 

oo 

o 

o* 

o- 

oo 

00 

00 

ao 

p-. 

r- 

o- 

0- 

00 

00 

OO 

CO 

P-. 

r-. 

4 — 

r-. 

ao 

OO 

CO 

CO 

P" 

r-. 

p-' 

r- 

00 

00 

co 

00 

r- 

r- 

r-- 

r-> 

-3 

w 

Q 

O 

x 

z 

o 

♦-4 

04 

04 

04 

04 

04 

CM 

04 

CM 

04 

CM 

CM 

CM 

04 

CM 

04 

CM 

CM 

CM 

CM 

04 

CM 

CM 

CM 

CM 

CM 

04 

04 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

04 

04 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

o 

O 

O 

O 

O 

O 

O 

O 

o 

o 

O 

O 

O 

O 

o 

o 

o 

O 

o 

O 

O 

O 

o 

0 

0 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

0 

0 

O 

O 

O 

O 

< 

uo 

O 

lA 

t-4 

lA 

o 

lA 

f— 4 

«A 

o 

lA 

r—4 

lA 

O 

lA 

r—4 

UO 

O 

lA 

r—4 

lA 

o 

UO 

r—4 

UO 

O 

UO 

r—4 

UO 

O 

UO 

H 

UO 

O 

UO 

r—4 

UO 

O 

UO 

r—4 

UO 

O 

UO 

r—4 

UO 

O 

UO 

r—4 

H 

►-4 

04 

04 

r— 4 

04 

r—4 

04 

r— 4 

CM 

*— 4 

CM 

r—4 

CM 

r — 4 

CM 

r—4 

04 

r — 4 

CM 

r—4 

CM 

r—4 

CM 

r—4 

CM 

r—4 

►—4 

o 

w 

UO 

UO 

UO 

LA 

UO 

lA 

lA 

lA 

lA 

UO 

lA 

UO 

xA 

lA 

lA 

lA 

o 

o 

O 

o 

O 

o 

O 

O 

O 

O 

0 

O 

O 

O 

0 

O 

0 

O 

0 

O 

O 

O 

0 

O 

O 

O 

O 

O 

O 

O 

0 

O 

a: 

PO 

m 

rO 

ro 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

co 

CO 

CO 

o 

o 

O 

o 

o 

o 

o 

O 

O 

O 

0 

O 

O 

O 

0 

O 

UO 

UO 

UO 

UO 

UO 

UO 

UO 

UO 

UO 

UO 

UO 

uo 

UO 

UO 

UO 

UO 

Ox 

ox 

O' 

Ox 

OX 

Ox 

Ox 

Ox 

Ox 

Ox 

Ox 

Ox 

Ox 

Ox 

OX 

Ox 

r-*- 

r- 

r- 

r- 

o- 

r>. 

p^ 

P^ 

P-* 

r- 

p"» 

r-' 

P^ 

r"- 

p"- 

O* 

lA 

UO 

UO 

UO 

UO 

UO 

UO 

UO 

UO 

UO 

UO 

UO 

UO 

UO 

UO 

UO 

• 

f— 4 

H 

H 

H 

H 

H 

r—4 

H 

«— 4 

f— 4 

r—4 

»— 4 

r— 4 

r—4 

r—4 

r—4 

co 

co 

co 

co 

CO 

CO 

CO 

CO 

CO 

co 

CO 

co 

CO 

co 

CO 

CO 

00 

00 

00 

00 

00 

00 

CO 

00 

00 

CO 

00 

00 

00 

00 

ao 

00 

w 

CO 

< 

H 

r—4 

r-4 

f— 4 

H 

r-4 

r— 4 

r-4 

r—4 

r—4 

f— 4 

•—4 

r—4 

r—4 

r—4 

r—4 

r—4 

H 

r— 4 

r—4 

r—4 

r—4 

r—4 

r—4 

r—4 

H 

r—4 

r—4 

r—4 

r—4 

•— 4 

r—4 

rH 

r—4 

r-4 

r—4 

r—4 

r — 4 

r—4 

r-H 

r—4 

r—4 

r-H 

r—4 

r4 

r-4 

r—4 

r—4 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

— 

X 

X 

X 

X 

X 

X 

o 

pH 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

O 4 

04 

04 

04 

CL, 

Or 

a* 

Om 

Pm 

Pm 

Pm 

Pm 

Pm 

P4 

Pm 

H  04 
CO 

Ui  Q 

f — 4 

H 

f — 4 

r-4 

H 

H 

r4 

r—4 

CM 

CM 

04 

CM 

04 

CM 

CM 

CM 

r—i 

r—4 

r—4 

r — 4 

r—4 

r—4 

r—4 

r—4 

CM 

CM 

CM 

CM 

04 

CM 

04 

CM 

r—4 

t— 4 

r—4 

H 

r-^ 

r-4 

r — 4 

r—4 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

H 

f- 

H 

H 

H 

H 

H 

H 

H 

H 

H 

H 

H 

h* 

H 

M 

H 

H 

H 

H 

H 

H 

H 

H 

H 

H 

f- 

H 

H 

H 

H 

H 

H 

H 

H 

H 

H 

H 

H 

H 

H 

H 

H 

H 

H 

H 

H 

H 

H  <5 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

to 

CO 

CO 

CO 

CO 

co 

co 

co 

CO 

CO 

CO 

CO 

co 

CO 

CO 

CO 

CO 

CO 

CO 

co 

co 

co 

co 

co 

co 

co 

co 

CO 

CO 

CO 

00 

CO 

CO 

CO 

to 

05 

Q 

Q 

Q 

a 

Q 

a 

a 

Q 

Q 

Q 

a 

o 

Q 

Q 

a 

Q 

a 

o 

Q 

Q 

Q 

Q 

a 

Q 

Q 

Q 

Q 

0 

Q 

Q 

a 

Q 

Q 

a 

Q 

Q 

Q 

O 

O 

Q 

a 

a 

a 

a 

0 

a 

O 

a 

35 


DST-1 


FREQ  RAINRATE  TFHP  PHASE  DIST  ERR  ABSORPTION  EXTINCTION  FIRST  EIGHT  LEGENDRE  COEFFICIENTS 

(GHZ)  (MM/HR)  DEG  K  COEFFICIENT  COEFFICIENT  123456 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOCM 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOoO 
oo  oooooooooooooooooooooooooooooooo<r 
ooooooooooooooooooooooooooooooooo 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOro 

ooooooooooooooooooooooooooooooooa> 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO' 

ooooooooooooooooooooooooooooooooo 


oooooooo 

oooooooo 

oooooooo 

oooooooo 


<Nr-400CN—4©0 

OOOOOOOO 

OOOOOOOO 

OOOOOOOO 


HO'(7'NfOOa'OJ4 
inrHOO^lCMOOr-* 
C  OOOOOOOO 
OOOOOOOOO 


H  1—4  *—i  -sj  r—i  r—i  r—i  kQ 

OOOOOOOOO 

OOOOOOOO 


oonr-fsirvcnNCJ 

CNr-400CN»-400 

OOOOOOOO 

OOOOOOOO 


QOmcoNajinooNoo 

CO»-400fnr*-400QO 

OOOOOOOOCN 

OOOOOOOOO 


^mmoor-'r-inco 

^CQNONvJCQ^OO 

r-400CNr-400m 

OOOOOOOO 


M  (N  CO  n  N  CN  <f 
(NHO^Mr^O^ 
OOOOOOOfO 


CNOOOOCOCNOOONCO 

Ovao<Nroo\oo<Nco 

MCNCNHrtCNOJH 

OOOOOOOO 


vO  ooonvoooonon 

co^^noOsj'anN 

OOOOOOOOO 


IflvOHCMiAOMOin 

00>OincO(Nh(S 

oovomor--'Ocom 

000*-<000.-4 


r**»oor^.mr^.oor^vo 

4coocN<r«oco 

Or'-'-tfiAOr-.-sfu-) 

«-4©0r-4r-400‘n 


h  in  in  O'  n 
HNHHiO  O  ^ 
CO  <0  O'  O  CO  on  ON 

•3  <3  ^  m  <3 


n  co  O'  m  in  03 

<N  r-'-  O  co 

O  n  co  O'  o 

in  <j  <3  ^3  in 


Oi  n  (N  O  H  O'  o 
»}  o  o  o  O'  fO  m 

CO  O'  O  N  o  ®  O' 

•3  *4-  m  -3  -3  *3  <3 


n  O'  on  in  0"0 
O'  H  rt  do  O  4 

00  O'  ON  O'  03  00 

3  3  3  <4-  -4-  3 


n  n  on  ®  in  cni  r-i 

n  3  O  3  r>.  <y  in 

CO  O'  CO  «  CO  O'  O' 

4  -o  *o  -O  4  oo 


Hvo>oo3voOfnOoJO'ri(nNO'HfO«0'NO*0(noO(on^OMn«jON'j 

ncgor-iw(»<finooO'00'cDO*00'00'40Nco*ONH>o^csfn«j^Nno' 

O'nofnH'OvfHincNjcocNjincNOJN-O'O'JHHiOcncocNiO'tninMO'cninco 

T-4rHOO*-4O©Or-4r-4OOr-4r-IOO*-4OOOCNr-4r-«OC0f-lr-4Oc0.-l.-<O0' 


«  l  •  l  i  •  i  (  t  i  I  i  i  <  i  I 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


ONNOrsNOininvon 
0O  H  W  N  H  O  in  t— 4 
NO'oocNJincn^m 
oinomr-coom 

'OO'O'HCMO'O'H 

in^00ini-400 


>onvoc44fnvocNr-( 
ncoo'OinpiO'^O' 
cn«Hrt(n  KHfOri 
OJCnJOO<J(NCnJCO*000 
«0O'C0H^O'00Hn 
m^ooin^oon- 


vooor^oocMCMinvo 

roO'OONpf^^in 
o\r*Hino<fcsr^r-i 
ooO'<fo>  o  <o  o\  co 
omiscor^mvoio 
r-»coor^r-co03 


vON'O'O'ON'ON 

3\ocN»n3vocNoo 

nOCO'OHnOOO'OM 

40Ninn<joiino 

OOnOvOOOnOO 

OfOO'O'OMOin 


on 

\D 

o 

m 

m 

VD 

CN 

co 

o 

vo 

o 

00 

O 

NO 

O 

00 

CO 

m 

m 

oo 

rv 

CO 

in 

NO 

r^- 

co 

oo 

o 

r- 

00 

00 

o 

3 

m 

3 

3 

r— 4 

co 

OO 

<o 

O 

3 

ON 

O 

>3 

ON 

n*  O 

00 

NO 

o 

in 

O 

oo 

co 

vO 

*3 

CN 

>3 

*3 

>3 

ON 

M3 

*3 

CN 

NO 

o 

o 

ro 

CM 

l-N 

r-* 

vo 

ON 

CN 

o 

co 

ON 

CN 

O 

CO 

VO 

r-4 

00 

vo 

CN 

co 

m 

CO 

00 

>3 

CN 

OO 

■3 

CN 

in 

3 

ON 

fN 

3 

m 

NO 

vo 

3 

CO 

co 

o 

r> 

QO 

CO 

o 

CO 

«— 4 

vO 

o 

CM 

3 

NO 

o 

3 

m 

r-4 

co 

00 

in 

r— l 

co 

00 

3 

m 

r-. 

«— 4 

O 

m 

rs. 

*— 4 

vo 

*>3 

r-4 

VO 

■*3 

r- 

H 

r— 4 

CO 

CO 

m 

NO 

in 

•3 

in 

m 

ON 

•3 

»n 

ON 

>3 

O 

<3- 

r-4 

o 

o 

-31 

*— 4 

O 

o 

co 

*— 4 

O 

O 

CO 

r-4 

o 

o 

o 

-3 

CN 

O 

O 

•3 

CM 

o 

r* 

CO 

r-4 

O 

rv 

co 

f— 4 

O 

co 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


HHHr4HHHH(NtMMCNN(NCSNHrlHHHrJHHNCJ(MCNCNMCSWH 


OOOOOOOO 


oooooooo 

cooooooor^r>.r>.r^ 

CMCN<N0M<N<NO4CN 


oooooooo 


oooooooo 

cooococooooo 

<N<NCM<N04<NCM<N 


OOOOOOOOO 


OOOOOOOOO 

flooococorsooow 

04CN<N04<NCN<N<NCN 


OOOOOOOO 


OOOOOOOO 

COCOOONONOOO 

<NCN<N<N<NCN<N<N 


OOOOOOOO 

OOOOOOOO 


inomr-Mmomf-i 

M  H  <N  H 


oooooooo 

oooooooo 


momr-4inoinr-4 

W  H  CN  r-4 


a oooooooo 

OOOOOOOOO 


inomi-iinoinHin 

<N  r-4  (N  H  CN 


oooooooo 

oooooooo 


OmHinoinnin 

r-4  CN  r-i  CM 


inininininminininininininmininooooooooooooooooo 

cococomrococococococorococococooooooooooooooooOin 


Xl 

nj 


0'0'(?'C^^O'0'0'0'^(7'0'0'0'C'(7'lN'rs-n«NNNrNNNh.Nr‘N(sNfvin 

HHHHr-tHHHHHr^HrMiHHHfnnncnc'C'fnfncntnnnnnfnfnco 


36 


85.50  10.00  280.0  1  10  .927913  1.716124  1.0000  .8244  .7779  .3927  .2139  .0903  .0388  .0147 


orvooH,-irNoa«QOHOa\ooH 

OOOOOOOOOOOOOO 


or'-aNvOOoor^vO>rMtr^iO'4<t 

O'cncnoOcnp-CNOvor-'CMfHvo 

^OO'J’CsJOvOOvjo^OO'vfO 

OOOOOOCMOOOCNOOO 


<jcomO<?'^r^r^r^ojr^r^f^.c*j 
N«jotn4invoinvOiovom\ovo 
mHO'CMftrlfOff'ONCIO'ON 
00O000>4»-“<«-M0<f»-<  — <  O 


Or^COiOnvOiOO^O'inO|A<7' 

vOOU3tnfO(N^sJOO'H<tOO' 

^\ONN«nvocoo''0(jieoo'\Da' 

oomcMOOvomcMOvocncMO 


ooor^r^r^mmcNcoomcNooc^ 

O(NOHCM(*Un^\0Mnr0\0N 

OvoocJCNMn^Nnin^Nfo 

cnHvOvJ^H^'XJvfOJO'vD^N 


O'JO'f'-'OHWvOHCT'CSVDHON 

oovvovoomcMr^OincMr-om 

0<JinCS<fr^vOCNlCOHvOCSCOH 

rsinavcofsinoJOcovoiNOcovo 


^0'NnOv0rl'fl0\NH\00'N 

oocNrsoOsjirtvocMOvJvooJO'J 

HWNHHHO'CT"OHO'0>\OH 

r-mooNOOvocMOONr^CMOOvr^ 


oooooooooooooo 

oooooooooooooo 

oooooooooooooo 

oooooooooooooo 


vO 

r- 

O' 

ON 

cn 

CM 

r- 

ON 

C-* 

CM 

r^- 

ON 

cn 

O 

vO 

o 

vO 

ON 

VO 

r- 

CM 

ON 

VO 

CM 

r- 

vo 

«0 

<n 

r-. 

CM 

vf 

m 

<r  cm 

-4 

m 

Mf 

O 

r- 

o 

4  CKO 

oo 

CM 

00 

CM 

Mf 

o 

00 

CO 

d 

00 

r>. 

r- 

CM 

00 

r—i 

CM 

oo 

cM 

o 

CM 

VO 

Cv 

CM 

t— < 

H 

vo 

CM 

H 

»— 1 

VO 

CM 

o 

CO 

H 

CM 

f-M 

CM 

O 

MT 

O 

r-» 

n- 

m 

m 

in 

CM 

ON 

r-~ 

m 

CM 

On 

OO 

vO 

CM 

m 

Mf 

r- 

Mt 

r- 

m 

o 

a\ 

VO 

oo 

ON 

CO 

vO 

>4 

r- 

00 

vo 

H 

un 

CM 

o 

in 

o 

m 

VO 

O 

m 

VO 

vO 

vO 

m 

r*. 

H 

o 

r-- 

CM 

H 

o 

CM 

in 

O 

00 

On 

m 

H 

H 

VO 

cn 

rM 

*— < 

vo 

cn 

H 

oooooooooooooo 


HHHHHHC'JCNOJNCNCNNC'J 


oooooooooooooo 


oooooooooooooo 

00Q0PvNNN«3O3«O0P^NNfs 

CNCNCMCMCMCMCMCMCMCMCMCMCMCM 


oooooooooooooo 

oooooooooooooo 


lAHiAOinHino/iHinomt-i 

CM  r-4  CN  O  NH 


oooooooooooooo 

^inmifli/itnmtAtninintninin 


lAiniAinviinintninininininin 

cocooooocooooococoooaooooooo 


37 


Table  8.  Input  for  Precipitation  Model  Test  Case 
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HEIGHT  PRESSURE  H20  DENSITY  ATTENUATION  COEFFICIENT  (DB/KM) 

CONDENSED  WATER 
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TOTAL  WATER  VAPOR  CONTENT  -  3.969  PRECIPITABLE  CM 
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ATMOSPHERIC  PROFILE  -  U.S.  STANDARD  (1962)  MODEL 
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GROUND  -  6  TEMPI  -  288.00  TEMP2  -  288.00  LP  -  .10000  LZ  -  1.00000  MV  -  .10000  E2RAN  -  (  .00000,  .00000) 

DEPRAN  -  50.000  H  -  .00000  DMV  -  .20000 

WT  -  .25000  GA  -  .45000  PO  -  .50000  ROCK  -  (  5.00000,  .10000) 


►H 

O 

u 

OS 

CO 

Cu 

O  Z 

o 

CM 

O  HM 

o* 

rH 

O  CO  0£$ 

VO 

VO 

O  CO  O 

VO 

• 

a  H 

r*** 

H 

X  U 

Ov 

CO  < 

Z  h 

H  cm 

cd  in 

^  hi  on 

> 

CO 

U  &  >~4 

w  p  CO 

ftS  CM 

o 

Cd 

Q 


O 

CM  PL, 

O 

H 

CM 

X 

H 

CO 

CO 

td 

X  z 

CM 

CM 

o 

O  H 

ON 

CM 

X 

r— 

H 

w 

O'  o 

CM 

VO 

3 

co  ►h  z 

H 

a$  5 

r* 

cti  co  O 

• 

vO 

H 

id  o 

m 

CM 

H 

O  Z 

O 

O  M 

O  CO  0<Z 

O  M  O 

M  H 

X  o 
co  < 
Z 


% 

/s 

O  • 

fd 

D 

h4 

O 

Z 

o 

< 

o 

^  Q 

O 

< 

<  o 

u 

in  cm 

Id 

H 

H 

Id 

in 

otS 

*~4 

Q 

Id 

e 

O 

H 

O 

H 

itf  Q 

55 


FREQUENCY  NADIR  ANGLE  EMISSIVITY  .  SURFACE  PARAMETERS 

(GHZ)  (DEG)  HORIZONTAL  VERTICAL  E2  EM 
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TEST  CASE  3A  -  VEGETATION  (MV-0.5)  -  TROPICAL  -  LIGHT  RAIN  -  NO  MS 
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ATMOSPHERIC  PROFILE  -  TROPICAL  MODEL 
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GROUND  -  5  TEMPI  -  300.00  TEMP2  -  298.00  LP  -  .10000  LZ  -  1.00000  MV  -  .50000  E2RAN  -  (  .00000,  .00000) 

DEPRAN  -  200.000  H  -  .10000  DMV  -  .20000 

W  -  .00330  VB  -  .10000  VBW  -  .05000  BULK  -  (  3.000000,  .000000) 

1  TEST  CASE  3B  -  VEGETATION  (MV-0.5)  -  TROPICAL  -  LIGHT  RAIN  -  MS 
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1  FREQUENCY  TOTAL  ATTENUATION  FOR  EACH  SPECIES 

(GHZ)  OXYGEN  H20  VAPOR  CLOUD  RAIN  TOTAL 
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6.  SUMMARY  AND  RECOMMENDATIONS 


6 . 1  Summary 

This  report  summarizes  enhancements  to  the  RADTRAN  transmittance/ 
brightness  temperature  computer  code  including  the  capabilities  to: 

(a)  evaluate  frequency  dependent,  polarized  surface  emissivity  -  a  menu 
driven,  user  selected  surface  type  solution  has  been  implemented  based  on  the 
surface  emissivity  submodels  described  in  Section  4.2.  The  surface  emissivity 
supports  evaluation  of  surface  emitted  brightness  temperatures;  (b)  calculate 
scattering  properties  of  precipitation  -  The  subprogram  for  calculating 
precipitation  optical  properties  described  in  Section  4.3  has  been 
incorporated.  This  provides  the  user  with  frequency  dependent  values  of 
extinction  coefficient,  single  scattering  albedo,  asymmetry  factor,  and  the 
angular  scattering  function  for  liquid  and  glaciated  precipitation.  These 
quantities  are  essential  to  performing  the  scalar  multiple  scattering  cal¬ 
culation  of  brightness  temperatures  in  the  presence  of  precipitation; 

(c)  perform  multiple  scattering  brightness  temperature  calculations  in  pre¬ 
cipitating  conditions:  an  exact  multiple  scattering  approach  for  fully 
polarized  brightness  temperature  calculations  has  been  included  as  a  user 
selectable  option  described  in  Section  4.4.  This  option  is  applicable  to  the 
simulation  of  data  from  sensors  with  polarization  discrimination;  (d)  perform 
statistical  retrievals  of  relevant  meteorological  parameters  based  on  RADTRAN 
simulations;  and  (e)  assimilate  user  specified  input  profiles  -  standard 
meteorological  profile  models  (e.g.  U.S.  standard,  tropical,  etc.)  has  been 
augmented  by  the  capability  to  accept  field  data,  radiosonde  and  other  upper 
air  data,  and  user  specified  format  profile  data  as  input. 

With  these  enhancements  and  the  algorithm  upgrade  described  in 
Section  4.1,  the  AFGL  RADTRAN  transmittance/brightness  temperature  computer 
code  provides  a  user  friendly  system  simulation  and  retrieval  tool  consistent 
in  purpose  and  format  with  other  AFGL  atmospheric  transmittance  models,  but 
suited  for  those  users  with  a  special  interest  in  the  microwave/millimeter 
wave  spectral  region. 
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6.2 


Our  major  recommendation  concerns  the  current  RADTRAN  gas  absorption 
formulation.  Comparison  with  FASCODE,  for  example,  does  indicate  a  difference 
in  the  calculation  of  gas  absorption  coefficient.  At  the  present  time, 
however,  it  was  decided  not  to  modify  the  RADTRAN  gas  absorption  models.  This 
decision  was  based  on  a  number  of  considerations  including:  (a)  maintaining 
model  consistency,  (b)  simplicity  of  the  current  approach,  and  (c)  technical 
direction  from  AFGL. 

In  light  of  the  requirements  for  accurate  gas  absorption  calculations 
imposed  by  physical  retrieval  procedures,  we  recommend  that  future 
consideration  should  be  given  to  replacing  the  current  RADTRAN  gas  absorption 
approach  with  one  similar  to  that  employed  in  the  FASCODE  "LASER"  option. 

This  would  employ  a  specialized  microwave/millimeter  wave  absorption  line 
parameter  file  for  the  region  such  as  that  used  in  the  FASCODE  "TEST”  option 
and  an  improved  layering  approach.  In  addition  to  improving  accuracy  and 
efficiency  of  the  model,  the  layering  update  would  provide  the  additional 
benefit  of  making  the  multiple  scattering  option  more  efficient  by  reducing 
the  number  of  layers  treated  within  the  precipitating  cloud. 

Finally,  for  the  purpose  of  physical  retrievals,  a  line-by-line  approach 
such  as  RADTRAN  is  much  too  inefficient  computationally  to  employ  for 
extensive  simulation  studies.  The  line-by-line  model,  however,  can  be  applied 
to  the  evaluation  of  appropriate  "band  model"  coefficients  to  be  applied  for 
specific  sensor  channels.  One  approach  for  the  microwave  has  recently  been 
suggested  by  Eyre  and  Woolf  (1988).  It  would  be  extremely  useful  to  provide 
this  code  along  with  the  RADTRAN  line-by-line  algorithm. 
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